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We investigate an interacting two-fluid cosmological model and introduce a scalar field representa- 
tion by means of a linear combination of the individual energy densities. Applying the integrability 
condition to the scalar field equation we show that this "exotic quintessence" is driven by an expo- 
nential potential and the two-fiuid mixture can be considered as a model of three components. These 
components are associated with baryonic matter, dark matter and dark energy respectively. We use 
the Simon, Verde & Jimenez{2005) determination of the redshift dependence of the Hubble param- 
eter to constrain the current density parameters of this model. With the best fit density parameters 
we obtain the transition redshift between non accelerated and accelerated regimes Zacc = 0.66 and 
the time elapsed since the initial singularity to = 19.8Gyr. We study the perturbation evolution of 
this model and find that the energy density perturbation decreases with the cosmological time. 
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I. INTRODUCTION 

Astrophysical data suggests that the Universe is ac- 
celerating [Ij, |3] . This acceleration may be explained by 
very different models, among them, the simplest one is 
the ACDM 0, It assumes a cosmological constant 
arising from the energy density of the zero point fluctu- 
ations of the quantum vacuum and cold dark matter in 
form of pressureless dust. While it fits rather well all 
the observational constraints, the small positive value of 
the energy density of the vacuum remains as an explana- 
tory challenge for physics today. See also Refs. [1, 111 for 
objections to this interpretation. The next step is to pro- 
pose a dark energy component that may vary with time 
and that is generally modeled by a scalar field. Most of 
these models assume that dark matter and scalar field 
components evolve independently. Again, this is not the 
solution because in the analysis of these models through 
SNIa or WMAP data, best fit models for one set data 
alone is usually ruled out by the other set at a large con- 
fidence limit i7|. As these conclusions are valid only for 
standard models, where dark energy and dark matter are 
decoupled, many papers have been devoted to interacting 
models [Sj. Certain models conceive the interaction as a 
time- variable dark mass, evolving with an inverse power 
law potential or an exponential potential, [Hj, 

In this paper we clarify this point and establish an 
exactly solvable model with a smooth transition from a 
matter dominated phase to a period of accelerated expan- 
sion. We introduce an interacting two-fluid cosmological 
model and investigate the effects of imposing the integra- 
bility condition of the whole equation of conservation. 
It forces the dynamic of the model to be governed by 
a modified Friedmann equation with three components. 
One of them, associated with an exponential potential, 
drives an exotic scalar field, (exotic quintessence). Here, 
we show that the problem of an accelerating universe can 



be realized in a comparatively simple manner within the 
framework of general relativity. Finally, the perturbation 
evolution of the model is investigated. 

Our paper is organized as follows. In section II we con- 
sider the general interacting two-fluid cosmological model 
and introduce the exotic quintessence. There, we obtain 
the evolution equation of the exotic field, find their im- 
plicit solutions, build the modified Friedmann equation 
and show the asymptotic behavior of the scale factor by 
using stability analysis. In section III we introduce bary- 
onic matter, dark matter and dark energy components 
and show that dark components satisfy separately an ef- 
fective equation of conservation with variable equation 
of state. In section IV we find confidence regions for the 
parameters of the model by using the Hubble function 
H(z) data, the age of universe and the redshift of the 
transition from non-accelerated to accelerated regime for 
the best fit model. In section V we present the equations 
governing the perturbations of the model. In section VI 
we express the conclusions. 



II. EXOTIC QUINTESSENCE 

Cosmological models are based in the Einstein equa- 
tions of the gravitational field where the source in- 
cludes different kinds of matter known, for instance pro- 
tons, neutrons, photons, neutrinos, etc., as well as non- 
relativistic non-baryonic cold dark matter and dark en- 
ergy. We embark in a less ambitious project by con- 
sidering a model consisting of two perfect fluids with 



an energy momentum-tensor Tik — T^^^ ~\- T^j^' . Here 

T^"^ = [Pn +Pn)uiUk +Pn9ik, whcrc pn and p„ are the 
energy density and the equilibrium pressure of fluid n 
and is the four- velocity. Assuming that the two fluids 
interact between them in a spatially flat homogeneous 
and isotropic Friedmann-Robertson- Walker (FRW) cos- 
mological model, the Einstein's equations reduce to two 
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algebraically independent equations: 

iH^ = Pi+P2, (1) 

Pl + P2 + 3iJ((l + Wi)pi + (1 + W2)P2) = 0, (2) 

where a(t) is the FRW scale factor and H{t) = a/a is 
the Hubble expansion rate. We introduce an equation 
of state for each fluid component Wn ~ Pn/ Pn n = 1, 2 
and for simplicity we assume that wi , W2 are constants 
and pi, P2 > 0. This simplified model leads to a reduction 
of the number of fundamental parameters required to de- 
scribe observations. It can be considered as an advantage 
from the computational point of view. We choose units 
such as the gravitational constant is set to SttG = 1 and 
c= 1. 

The whole equation of conservation ([2]) shows the inter- 
action between both fluid components allowing the mu- 
tual exchange of energy and momentum, meaning that, 
there will be no local energy-momentum conservation 
for these fluids separately. Then, we assume an over- 
all perfect fluid description with an effective equation of 
state, w = p/ p = —2H/3H^ — 1, where p = pi + P2 and 
p = pi + P2- So that, from Eqs. fl])-© we get 

- 2ij = (1 + Wi)pi + (1 + W2)P2 = (1 + w)p. (3) 

To avoid an eventual super acceleration of the universe, 
which could lead to a "big rip" singularity, we choose 
wi > -1 and W2 > -1 [11. Hence, {l + w)p = -2H > 
and there is no big rip singularity. These models can be 
investigated by introducing a scalar field representation 
of the interacting two-fluid mixture 



0^ = (1 + Wi)pi + (1 + W2)P2, 



(4) 



with (j)'^ = —2H. The dynamic equation for the scalar 
field is obtained from the equation of conservation ([2]) 



-(l + Wi)H,j>- 



Wi - W2 P2 



= 0. 



(5) 



It can be integrated by setting the interaction between 
the two fluids by 



P2 + A(j)p2 = 0, 



(6) 



where ^ is a new constant parameter of the model. In- 
tegrating Eq. ([6]), we find that the energy density of 
the second fluid can be associated with an exponential 
potential 



P2 = P2aHQe 



(7) 



where p2o is a positive integration constant and Hq, (jy^ 
are the present values of Hubble constant and scalar flcld. 

From Eqs. ([5|) and (O we obtain the total en- 

ergy density and pressure of the fluid mixture and the 
dynamical equation for the scalar field 



Wi - W2 



Wl 



1 



V, 



(8) 



-{l + wi)Hcp- 



I Wl — W2 

Wl 1 + Wl 

Wl — W2 dV 



(9) 
(10) 



with p + p — 4>^ . These equations are different than the 
conventional ones describing quintessence, in contrast, 
define an exotic scalar field. When the interacting two- 
fiuids system is related to the scalar field in the form 
pi = 0^/2 and p2 = V{(f>), with equations of state pi = pi 
and p2 = —p2, meaning that wi — 1 (stiff matter) and 
u'2 = — 1 (vacuum energy)^ the exotic scalar field reduces 
to quintessence. Then, due to the interactions between 
the two-fluid components the energy-momentum tensor 
conservation of the system, as a whole, is equivalent to 
the Klein-Gordon equation. For any other interacting 
two-fluid mixture, the cosmological model contains an 
exotic quintessence field driven by an exponential po- 
tential. 

Using the integrability condition ([6]) in the field equa- 
tion (l5|), its first integral is given by 



and 



(j) = (j)^-A\n{l + z)-cHo [ 

Jo 



(1 + z)(l+3"'i)/2 



H 



(11) 



dz, (12) 



where c is an arbitrary integration constant, z = —1 + 
flo/a is the redshift parameter and oq is the present scale 
factor. This model is finally closed when the Eq. ([TT|) 
is inserted into the energy density ([5]) [l3| . Hence, the 
Friedmann equation ([1]) reads 



QcAHqH{1 + z 



,3(l+«)i)/2 



+ 



3(l + u;i) -^2 

3(wi - W2)P2 + ic^Hlil + z)3(i+-i)](13) 

As a consequence of the linear term in the expansion 
rate H , this equation can be seen as amodified Fried- 
mann equation. Its solution gives the scale factor and 
the model we propose, containing exotic quintessence, 
could be formally solved. 

In order to obtain the asymptotic behavior of the scale 
factor it will be useful to find the constant solutions of 
the dynamical equation for the overall equation of state 



w = -m{w-wi)i\j^{l + w)-{l + w)\, (14) 



and investigate their asymptotic stability. This equa- 
tion has two stationary solutions wi and A^ /"i — 1 (Eq. 
([3]) excludes the solution w — —1). Assuming that 
A^ < 3(1 -f Wl), we find that Wi is an unstable solu- 
tion while A^ /Z — 1 becomes asymptotically stable. Es- 
sentially, the evolution of the geometry is dictated by 



3 



w = —2H/3H^ — 1, meaning that the universe begins 
to evolve from an unstable phase as it were dominated 
by the first fluid wi at early times, a w ^2/3(i+mii)^ 
ends in a stable expanding phase dominated by the ex- 
ponential potential, a ~ t^^"^ . The latter becomes an 
expanding accelerated phase when the slope of the po- 
tential satisfy the inequality < 2. 



III. DYNAMIC OF BARYONIC AND DARK 
COMPONENTS 

The integrability condition ^ can be considered as 
an effective equation of conservation for the second fluid. 
This allows us to identify p2 with the energy density of 
the dark energy component and Wde = —1 + A(j)/3H with 
its effective equation of state. Expressing the latter in 
term of the exotic field (j) we get 



Wde 



-1 



A(f) 



-1 



A4) 



3 [02 _ W2V] 



and the relation 



-1 + —{1+Wdef, 



(15) 



(16) 



linking the overall and dark energy equations of state. 
Also, for convenience we write p2 = p20-ffo (1 + zY'''^ with 



3Aln(l + z) ^ -A{(t)- (t>o). 



(17) 



The H linear term in the Eq. is adequate to describe 
non-relativistic non-baryonic cold dark matter compo- 
nents whose energy-momentum tensor is approximately 
dust-like. Finally, the baryonic matter is introduced by 
setting wi = in the third term of the Eq. (fT5|) . Mak- 
ing these identifications the modified Friedmann equation 
(fT3l) becomes 



^ = ^drn{l + zf'^^+nde{l + zf^ + nk{l + z)\ (18) 
tin tlQ 



^0 

where 



2cA 
3-A^' 
-W2P20 
3-A^' 



3-A^' 



(19) 
(20) 
(21) 



are the present dark matter, dark energy and baryonic 
matter density parameters respectively. As these are con- 
strained according to fldm + ^dm + — 1, we conclude 
that 



(c + A)^ - W2P20 = 3. 



(22) 



— 1 < i«2 < —1/3. Hence, one finds that the conditions 
A^ < 3 and > are consistent with the requirement 
of having positive density parameters. 

The exotic quintessence is essentially based on the in- 
tegrability condition ([6]) under which the conservation 
equation ^ and the exotic field equation ([5]) can be in- 
tegrated. Actually, from the above condition cA > 
along with Eq. (|lip . we get A(j> > for an expanding 
universe. Then, p2 is a Liapunov function and the solu- 
tion p2 = 0, of Eq. is asymptotically stable. Also, 
from Eqs. ([7]) and pT|) the energy density of the first 
fluid has a vanishing limit in the remote future. Hence, 
this general model is viable and it does not contradict 
basic cosmological conjectures. 

From Eqs. and (I^T]) we can express the parameter 
A and the integration constant c in terms of the present 
density parameters ^dm, 



A^ 



± V3 r^dm 

± 2V3 



besides. 



W2P20 = 



12nde^b 



^dm 



(23) 
(24) 

(25) 



Finally the original problem of the interacting two- fluid 
mixture governed by system equations (Il|)-(I1]) is equiva- 
lent to an effective model with a "three-fluid" mixture. 
So that, the effective dynamical equations of our model 
read 



3iJ =Pdm+Pde+Pb, 
im + 'SH{1 + Wdrn)Pdm = 0, 
Pde + 3-ff (1 + Wde)pde = 0, 
Pb + 3Hpb = 0, 



where 



Pdn 
Pde 



= SHondmil + zf^^H, 
3H'Ml + zf^'+'\ 

Pb = + z)\ 



(26) 
(27) 
(28) 
(29) 



(30) 
(31) 
(32) 



are the effective energy densities of dark and baryonic 
components and 



2nb 



Wde 



^dm 



-2+^^^(l + z)3/2 
H 



(33) 



is the effective equation of state of the dark energy. Also, 
we find the following relation 



Wdn 



-1 



A^ 



(l+Wde)^ 



(34) 



During the accelerated epoch of the universe a > 0, the 
SEC is violated and p+3p = Pi + {l + 3w2)p2 < leads to 



between the effective equations of state of dark compo- 
nents. So that, the knowledge of them determines the 
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remaining ones including w. To obtain the above effec- 
tive dynamical equations of the model we have taken into 
account that pde oc p2- Hence, we have identified the 
equation of conservation ^ with to express, after 
using ([8]), the effective Wde for dark energy in terms of 
the observed density parameters, the present Hubble ex- 
pansion rate and the redshift parameter. 

In our model the stationary solutions w'^ = at early 



times and w 



-1 + at late times along 



with Eqs. p6p and p4|l lead to the stationary solutions 



-1 + A/V^, 

12 , 



-1 + ^2/3 and 



0, 



'^dm ~ + ^'^/3)/2. As the observed density parame 
ter satisfies the condition fi^^^ < 817;,, then < 2 by us- 
ing Eq. (|23p . In this case, we find that w'^^ is an unstable 
solution at early times and becomes asymptotically 
stable at late times. Here, the evolution of the geome- 
try represents a universe that begins to evolve as it were 
matter dominated at early times and ends in an accel- 
erated phase dominated by the dark energy component. 
As the accelerated epoch begins ata^Oorui^— 1/3, 
the corresponding redshift is given by the expression 



^lde{V2-Af 



(35) 



On the other hand, the Eq. ([35| allow us to find the 
elapsed time to from the creation 



to 



V3 
A 



(36) 



where we have used that the scale factor behaves as a ; 
ao{t/toy^^ , at early times. 



IV. OBSERVATIONAL CONSTRAINTS 

In Ref. [l^ it was used the recently published Hub- 
ble function H{z) data [SV&J(2005)] [l^, extracted from 
differential ages of passively evolving galaxies. This is in- 
teresting for, among other reasons, the function is not in- 
tegrated over, in contrast to standard candle luminosity 
distances or standard ruler angular diameter distances 
Since the Hubble parameter depends on the differen- 
tial age of the Universe as a function of z in the form 
H{z) = —(1 -I- z)~^dz/dt, it can be directly measured 
through a determination of dz/dt. In the procedure of 
calculating the differential ages, Simon et al. have em- 
ployed the new released Gemini Deep Deep Survey fv?] 
and archival data [ISi], (19] to determine the 9 numeri- 
cal values of H{z) in the range < z < 1.8, and their 
errors [see Table 1]. These data will be inserted in our 
Eqs. (|17p and to derive restrictions on the range of 
possible values for the density parameters p^ - (PT|) . 

We adopt the prior Hq = 72kras~^Mpc~^ for H{z — 
0) . It is exactly the mean value of the results from the 
Hubble Space Telescope key project [lO] and consistent 
with the one from WMAP 3-year result [2l| . 



z 


H{z) 


1(7 






uncertainty 


0.09 


69 


±12 


0.17 


83 


±8.3 


0.27 


70 


±14 


0.40 


87 


±17.4 


0.88 


117 


±23.4 


1.30 


168 


±13.4 


1.43 


177 


±14.2 


1.53 


140 


±14 


1.75 


202 


±40.4 



TABLE I: SV&J Hubble Parameter vs. Redshift Data. 
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Fig. 1: Observational H{z) with la uncertainties 
from SV&J(2005) and the best fit theoretical H{z). 



The parameters of the model can be determined by 
minimizing the function 



9 

X^i^dmMde) = ^ 

1=1 



[Hthi^dm, ^de: Zi) — Hob{Zi 



(37) 



where 



Hth{^dm,^de;Zi) = Hoil + Zif^^ X 




dm 



^ ± [(1 ± z.)3(^-i) - 1] 



(38) 

is the predicted value for the Hubble parameter, obtained 
from Eq. (fT8|) . and A is calculated from Eqs. (fT2l) and 

m- 

Hob is the observed value of H at the redshift Zi, is 
the corresponding Icr uncertainty, and the summation is 
over the 9 observational H{zi) data points at redshift 
[12]. Also, we adopt the prior = 0.05 [H]. 

We find a local minimum of for = 0.345865 and 
nde = 0.606504 (x^ = 8.69263) and present the observa- 
tional H{z) data in Fig. 1 with error bars and the theo- 
retical line corresponding to the best fit parameter. The 
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Acceleration 




z 



Fig. 2: Acceleration vs. redshift for the best fit 
model. 




OPEN MODEL 



-0.2 [ ai„.= 0,346 0,1. = 0.606, 
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Odm 

Fig. 3: The la and 2a confidence regions, (inside and 
between the elliptic contours) for Qdm and Qde from 

SVSzJ{2005). The cross is the best fit model. The 
straight line corresponds to the flat cosmology making 

the separation between open and closed universes. 
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Fig. 4: State parameters for dark matter, dark energy 
and overall fluid. 
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Fig. 5: Energy densities (in units of 3Ji"o) vs. redshift 
z 
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Fig. 6: The age of the universe (in units of Hq^). 
The solid line corresponds to our model and the 
dashed line represents the ACDM one. 

Fig. 2 shows that the Universe begins to accelerate on 
z ~ 0.66. A similar result can be obtained from ACDM 
flat cosmology when the density parameters are f2m = 0.3 
and Q\ = 0.7 [23]. In Fig. 3 we plot confidence regions 
in the ridm - ^de plane. The true values of those param- 
eters are inside the inner eUipse or between both ellipses 
with 68.3 or 95.4 percent of probability respectively. In 
Fig. 4 we plot the equations of state for dark matter 
Wdm = Pdrn/Pdm, dark energy Wde = Pde/pde and overall 
fluid w = p/ p. According to the previous stability analy- 
sis predictions we see that the asymptotic behavior of the 
overall state parameter varies from w ^ (near cold mat- 
ter behavior) in the far past to w = — 1 -I- /3 ^ —0.63 
in the far future (dark energy behavior). In Fig. 5 we 
show the energy densities corresponding to the baryonic 
and dark components in terms of the redshift z. 

We calculate the age of the universe ([55)1 with the best 
flt parameters and find that to = 19.885Gt/r. In Fig. 
6 we plot the time elapsed (in units of Hq^), since the 
initial singularity to present days, for our model and the 
flat ACDM model, as a function of the matter density. 
We also show the border to = llGyr coming from the 
bound of the oldest stellar ages. The age of the universe 
in this coupled scenario tends to be much higher when 
compared with the ACDM case [l^ 



V. LINEAR PERTURBATIONS 

Cosmological models with two interacting fluids have 
been investigated with the purpose of describe the evolu- 
tion of dark components. There the energy momentum 
tensor of the interacting components is not separately 
conserved. Usually these cosmological model are pre- 
sented with interacting matter species have a non con- 
stant equation of state parameter [2^ or with DE having 
a constant equation of state parameter coupled to DM 
[131, [l^. However, in the model we are investigating the 
interaction between the two fluids is setting by Eq. 
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For this choice the field equation can be integrated, gen- 
erahzing the case of quintessence driven by the exponen- 
tial potential. Then, it will be interesting to investigate 
the evolution of the density perturbation. 

In the synchronous gauge the line clement is given by: 



(39) 



where the commoving coordinate are related to the 
proper time t and position r by dr = dt/a^ dx — dr/a 
and hij is the metric perturbation. The scalar mode of 
hij is described by the two fields h{k,T) and r]{k,T) in 
the Fourier space. 



hij{x,T) 



d^ke 



k-x 



(40) 

with k =^ kk. The Einstein equations to linear order in 
k-space, expressed in terms of h and rj, are given by the 
following four equations [29]: 



la' 



h' 



' 2 a 
k^rl = 47rGa2(p 



2—h' - 2k^T] = -SirGa^ST^', 



h" + 6rj" 
-24TTGa^{p + p)a. 



2—(h' + 67j')-2k^T]- 
a 



(41) 
(42) 
(43) 

(44) 



Here, the quantities 9 and a are defined as (p + p)0 = 
ik^ST^, ip + p)a ^ -{k,kj - 5,j/3)Sj and = T] ~ 
SjTj^ /3 denotes the traceless component of the tensor Tj. 
In addition, 9 is the divergence of the fluid velocity 6 = 
ik^Vj and ' means d/dr. 

Let us consider a fluid moving with a small coordinate 
velocity = dx^ /dr, then, can be treated as a pertur- 
bation of the same order as energy density, pressure and 
metric perturbations. Hence, to linear order in the per- 
turbations, the energy-momentum tensor, with vanishing 
anisotropic shear perturbation , is given by 



t; 



T' 



5p), 
{p + 6p)6}. 



j 



(45) 
(46) 
(47) 



For a fluid with equation of state p = wp, the per- 
turbed part of energy- momentum conservation equations 



T.^"^ = in the k-space leads to the equations 
6' ^-{l + w){9+'^)-3H{^^^-w)6, 



9' = -nil-3w)9- 



Sp/Sp 



1 



w 



1 



w 



k^5, 



(48) 
(49) 



where S = Sp/ p and Ti. ~ a' /a ~ aH ~ a. Besides, using 
equations (gl]), and (gH) we arrive at 



h" + Hh' + 3H^ [1 



(50) 



We have showed that our interacting two-fluid model 
can be associated with an overall perfect fluid descrip- 
tion based in an effective equation of state w = {wipi + 
W2P2)/{pi + P2)- Hence, we investigate the asymptotic 
regimes at early and late times assuming nearly constant 



w 



and w 



w 



-1 + AVa 



equations of state w ■ 
respectively. 

At early time, when the overall fluid has w ~ 0, the 
effective fluid perturbations evolve similar to those of or- 
dinary dust with 9 = 9 = 0, and from Eqs. (|48ll50p we 
obtain 



5 + 2H5- -H^S = 



(51) 



and 5 =^ cit^^ + C2t^^^, where ci and C2 are arbitrary 
integration constants. In this dust dominated era the 
perturbation grows as 5 « a showing an initial unsta- 
ble phase and compatible with the observation that the 
primordial universe would have tiny perturbations which 
seed the formation of structures in the later universe. 

At late times, we are interested to find the evolution 
of the linear scalar perturbations for any mode k. To 
this end we write the second order differential equation 
for the density perturbation 5 and the first order differ- 
ential equation for the divergence of the fluid velocity 9, 
evaluating them on the asymptotically stable equation of 
state w ~ wK In this case, from Eqs. (|48ll50p we get: 



6" + m' + 



w^k^ ~^{l + w^){l + iw^)re 



+3w'(l + w')H6l = 0, (52) 
^' = -W(l - 3w> + (53) 

Taking into account that in the late time regime the 
scale factor behaves as a oc t2/3(i+«; ) ^^^^ calculate 
the conformal time r, a and Ti. = a' j a 



^ ^(l-l-3«.')/3(l+u.') 



^2/(l+3i«') 



u = 



(1 + 3u;0t-' 



(54) 
(55) 

(56) 



From Eqs. (|52p and (|53p the perturbation evolution be- 
comes mode dependent with the fc^/7i^ term, and for low 
energy modes their solutions can be obtained assuming a 
power law dependence of the perturbations with the scale 
factor, (5 cx a" and 9 oc a^. In this case the approximate 
solutions for = —1 + /3 = —0.63 are given by 



(57) 
(58) 



j2.89 

Si 
^0.55 



^0.S9 



^3.33 ■ 



where 9q, 5i and 62 are integration constants while 9i is 
a function of 60 and wK This shows that the coupling 
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to 9 in Eq. ([52)) can be neglected for all scales we are 
interested. Finally, expressing the Eq. ([5^ in term of 
conformal time we get 



5" 



5' 



1 + Sw' T 



wk — 6 ■ 



1 



1 



1 + 3uf 



S^O. (59) 



The general solution of the latter equation in terms of 
the Bessel functions is 



with 



-l + 3u;' 
2(l + 3?«')' 



± 



2(1 + 3m;0' 



(60) 



(61) 



At late times, it can be approximated by the two first 
terms of the Eq. ([58| showing that the energy density 
perturbation decreases for large cosmological times for 
modes satisfying the condition k^/H^ <C 1. For high 
energy modes, k^/Ti,'^ 3> 1, the perturbation 



1 



(62) 



decreases but slowly that the low energy modes. This 
results can be understood writing the Eq. ([52| as the 
equation of motion of a dissipative mechanical system by 
using the analogy with the classical potential problem 



d7 



2 



where 



V{5) = w'k^ 1 



ntJ 2 



(63) 



(64) 



D{6, 5') = -{l + w^){l + 3w^)HH'5^ + H6'^, (65) 



n_/2 
/to 



2w'fc 



3(l + w')(l + 3w')' 



(66) 



The potential V has an extreme at S — 0, it is max- 
imum for Ti < Tio or a minimum for Ti. > Tig- On 
the other hand, assuming that the perturbation depends 
on the scalar factor in the form 6 cx a", we find that 
D w 0.087i^(5^ > 0. Hence, for any mode k the pertur- 
bation begins to grow at early times for Ti. < Hq, while 
at late times for Ti. > Tia, the function inside the square 
bracket in Eq. (j63p is a Liapunov function and the per- 
turbation decreases asymptotically reaching ^ = in the 
limit t — !■ CX). 



VI. CONCLUSIONS 

We have shown an interacting two-fluid cosmological 
model that allows us to reproduce the accelerated behav- 
ior of our universe and its probable age . The model gives 



rise to an exotic scalar field dubbed exotic quintessence 
which reduces to quintessence when one fluid is associ- 
ated with stiff matter and the other with vacuum energy. 
Setting the interaction between the two fluids by Eq. ©, 
the field equation is integrated, generalizing the case of 
quintessence driven by the exponential potential and, the 
equation governing the scale factor looks like a mod- 
ified Friedmann equation. 

We have obtained the evolution equation for the over- 
all equation of state of the model and showed the asymp- 
totic behavior of the scale factor i.e., the universe be- 
gins from an unstable phase dominated by the first fluid, 
a « ^2/3(i+«ii)^ ends in a stable expanding phase 
dominated by the exponential potential, a « t^^^ . The 
latter becomes accelerated when the exponential poten- 
tial slope satisfy < 2. Setting wi = 0, the scale factor 
interpolates between pressureless matter and dark energy 
phases. 

Using the Hubble function H{z) data from Table 1 
we minimize the function ()37p . obtaining the best 



fit densities parameters ildm = 
0.606t^:Ej26 with a reduced 



0.346 



+0.054 
-0.046 



and ilj 



1.24. These results are 
consistent with those found in the literature, see for in- 
stance Ref. [30] for null coupling {ridm — 0.37j^o qs) 
with the result obtained in Ref. [Slj, (Odm = 0.30^0 07) 
through mean relative peculiar velocity measurements for 
pairs of galaxies. With our best densities parameters and 
the priors for Hq — 72kms~^Mpc~^ and fib = 0.05, we 
obtain the theoretical H{z) function, plotted together 
with the SVJ(2005) experimental data in Fig. 1. In Fig. 
2 we plot the acceleration of the model as a function of the 
redshift and find the transition from the non accelerated 
phase to the accelerated one around Zacc = 0.66 '24]. Our 
value agrees with the result obtained by a nearly model 
independent characterization of dark energy properties as 
a function of redshift, {zacc = 0.42tJ] [56,[3). The prob- 
lem of why an accelerated expansion should occur now 
in the long history of the universe seems to be naturally 
dressed in our model. Considering the age of the uni- 
verse, we take into account that the age of the oldest stel- 
lar objects have been constrained for instance, by using 
a distance-independent method js^l, (^o = 13.5 ± 2Gyr 
for Globular clusters in the Milky Way) and the white 
dwarfs cooling sequence method [3J] {t = 12.7 ± O.lGyr 
for the globular cluster M4). Then, the age of universe 
needs to satisfy the lower bound to > 12 — 13Gyr. This 
condition is fulfilled by our model with to = 19.885G'yr, 
as it can be seen in Fig. 5. 

The energy density perturbation of the model grows 
in the first stage of the universe showing that initial in- 
stabilities in the primordial universe could leads to the 
formation of structure in the later universe. At late times 
we have found a Liapunov function which indicates that 
the perturbation decreases asymptotically reaching (5 = 
in the limit i — > 00. 
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